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Abstract Analysis of an exactly soluble model of phonons 
coupled to a carrier in a quantum dot provides a clear illus- 
tration of a phonon bottleneck to relaxation. The introduc- 
tion of three-phonon interactions leads to a broad window 
for relaxation by the processes of LO phonon scattering and 
decay. 

1 Introduction 

The possible existence of a phonon "bottleneck" that 
might inhibit the relaxation of carriers in quantum dots 
remains a subject of controversy. Although the bottle- 
neck is readily eliminated from theoretical treatments 
by the introduction of non-phonon scattering processes, 
there has also been interest in demonstrating relaxation 
by intrinsic phonon scattering mechanisms alone. 

Within the one-phonon approximation to the self- 
energy, the carrier's spectral function (density of states) 
consists of a number of peaks which broaden rapidly with 
increasing temperature. This level broadening is known 
not to be a lifetime effect [Q, but might be supposed to 
help relieve the constraint of exact level matching in LO 
phonon emission. Indeed, the non-vanishing relaxation 
rates found in simple non-equilibrium treatments ^,|| 
can be traced to the non-vanishing spectral widths of 
the one-particle Green's function. 

Though physically plausible, such theories conflict 
with the results of experiments by predicting a vanish- 
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Fig. 1 Electronic spectral function at 300 K for an unper- 
turbed level separation A = hu>o = 36meV. The electron- 
phonon interaction is responsible for the Rabi splitting of 
the peak for the first excited state. 



ingly small relaxation rate in the limit of low tempera- 
ture and low excitation density, where Auger processes 
can be neglected. For this reason we have made a one- 
phonon calculation of the carrier's spectral function in an 
exactly soluble model for which we know that there can 
be no relaxation, as will be explained in Sect. ^| below. 
Figure |l| illustrates the broadening of the spectral peaks 
that — in this case, at least — is a spurious result due to 
the one-phonon approximation. Accordingly, one should 
regard with caution any proposed relaxation mechanism 
that relies on the broadening obtained in this widely- 
used approximation. 

Despite the impossibility (in our model) of obtain- 
ing relaxation by LO-phonon emission alone, we show 
in Sect. |^ that the introduction of three-phonon inter- 
actions (leading, e.g., to the decay of an LO phonon to 
give two LA phonons) can give a broad, temperature- 
dependent window for relaxation, as was found previ- 
ously Q by use of the Wigner-Weisskopf method. 

2 A soluble model 

Our model for phonons interacting with an electron in 
a parabolic confining potential is the usual electron- 
phonon Hamiltonian with a special choice of matrix ele- 
ments \/n + 1 g q between neighbouring electron levels, 



H 



n q 

Vn + lg q (b q + bl q ) (4 



E 
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(1) 



where c n and b q are the electron and phonon annihilation 
operators. The calculation of the spectral function illus- 
trated in Fig. [jjwas based on Eq. (jl]), truncated to four 
levels. In the one-electron subspace, the ladder operators 
a = Y2 n \/n + 1 c^Cn+i and obey Bose commutation 
rules and the above Hamiltonian is equivalent to 



H = Aa^a - 
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{%u Jq b\b q + g q {b q + bl q ){a + a))}. (2) 



The Hamiltonian (||) represents a system of coupled os- 
cillators, so that all features of its energy spectrum and 
response can be obtained without approximation. For 
example, the response to an impulsive electric field ap- 
plied at time t — is given by the retarded susceptibility 



x(t) - -i(x(t)x(o) - x(p)x(t)) e(t) , 



(3) 



where X 



is proportional to the electron dipolc 



moment and 6(t) is the unit step function. In the fre- 
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quency domain, x satisfies the Dyson equation 

X ( W ) = x » + xV) J2 q \g q \ 2 D» xM , (4) 

which can be solved for x(w) in terms of the unperturbed 
phonon Green's functions D®(w). 

For coupling to dispersionless LO phonons the calcu- 
lation can be completed analytically. The spectral func- 
tion — 2Imx(w) consists of a pair of (Rabi split) delta- 
function peaks that result from the binding of a phonon 
to the electron jj]. In the time domain, x(t) is the sum of 
two undamped oscillations: the electric field excites the 
electron into a state which is a superposition of exact 
eigenstates from the discrete spectrum of (|J), and from 
which the electron is unable to relax. 

For coupling to LA phonons, the sum appearing in 
can be evaluated numerically. We find that if A falls 
in the LA phonon continuum, relaxation is possible by 
emission of an LA phonon, but the rate is many orders 
of magnitude too slow for single LA-phonon emission to 
be significant as a mechanism for relaxation. 

3 Influence of three-phonon interactions 

Finally we turn to mechanisms that can lead to the fast 
relaxation of an electron in a quantum dot. Interactions 
between phonons of branches i, j, and k can be described 
by a phenomenological Hamiltonian || 

H^ k =/3 J2 44"rK X i X r> ( 5 ) 
p+q+r=0 

where X* = (&* + b l J p ) / (M l to p ) V 2 . We choose the coef- 
ficient (3 (the same for all ijk) to obtain agreement with 
the observed lifetime of free LO phonons with respect 
to decay to two LA phonons, r = 2.5 ps in GaAs at 
300 K. Note that we must, more generally, consider also 
the interactions with i = LA and j = k = LO, because 
an LO phonon bound to an electron can be liberated by 
absorption or emission of an LA phonon with simultane- 
ous relaxation of the electron. In contrast, non-dispersive 
free LO phonons cannot be scattered in this way because 
of the constraint of energy conservation. 

The effect of the interactions is calculated perturba- 
tively. The LO-phonon Green's function D q (u>) satisfies 
a Dyson equation 

D q {u) = ti° q (cj) + D° q (cj) n q {u) D q (cj) , (6) 

where the phonon self-energy n q {uj) is evaluated to sec- 
ond order in /3 by use of the Matsubara technique ||. 
The phonon Green's function obtained in this way re- 
places Dq(u>) in Eq. (^), which gives an approximation 
for the susceptibility for the electron. A relaxation rate 
for the electron is estimated by fitting an exponential Q 
to the time-domain susceptibility, x(*)- 

Results of these calculations are shown in Fig. |[ for 
temperatures T = and 300 K, and for a range of elec- 
tron level separations. In all cases the relaxation rate is 
greatest when the unperturbed level separation matches 
the LO phonon energy, A ~ hug. This is understandable 
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Fig. 2 Carrier relaxation rate as a function of level separa- 
tion. Dashed lines are results for the mechanism LO — + 2 LA ; 
results shown by solid lines also include LO — > LO ± LA. 

on qualitative grounds, as in this region the electron- 
phonon bound-state wave function contains nearly equal 
amplitudes for the electron to be in its ground state 
(with one additional phonon) as in its first excited state, 
so that the interaction Hamiltonian has a large matrix el- 
ement between the bound state and the electron ground 
state. A new result of our work is that the relaxation rate 
remains high only up to a threshold value of A beyond 
which the energy conservation law cannot be satisfied; in 
fact, as can be seen from the graph, there are slightly dif- 
ferent thresholds for the two kinds of relaxation process 
considered. 

The scattering process LO — > LO — LA has a partic- 
ularly strong dependence on temperature. Below reso- 
nance, it leads to no enhancement to the relaxation rate 
at low temperature because there are few LA phonons 
available to be absorbed. In contrast, at 300 K the ab- 
sorption of LA phonons enhances the relaxation rate by 
more than an order of magnitude at the smallest values 
of A, and so provides the major route for relaxation in 
this regime. 
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